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ABSTRACT 

We present numerical simulation of bulk Lorentz factor of relativistic electron-positron 
jet driven by Compton rocket effect from accretion disc radiation. The plasma is as- 
sumed to have a power-law distribution n e (pi) oc 7 _s whith 1 < 7 < j max and is con- 
tinuously reheated to compensate for radiation losses. We include full Klein-Nishina 
(hereafter KN) cross section, and study the role of energy upper cut-off "f m ax, spectral 
index s, and source compactness. We determine terminal bulk Lorentz factor in the 
case of supermassive black holes relevant to AGN and stellar black holes relevant to 
galactic microquasars. In the latter case, Klein-Nishina cross section effect are more 
important, and induce terminal bulk Lorentz factor smaller than in the former case. 
Our result are in good agreement with bulk Lorentz factors observed in galactic sources 
(GRS1915+105, GROJ1655-40) and extragalactic ones. Differences in scattered radi- 
ation and acceleration mechanism efficiency in AGN environment can be responsible 
for the variety of relativistic motion in those objects. We also take into account the 
influence of the size of the accretion disc; if the external radius is small enough, the 
bulk Lorentz factor can be as high as 60. 

Key words: galaxies: active - galaxies: jets - radiation mechanism: miscellaneous - 
stars: individual: GRS1915+105, GROJ1655-40 



INTRODUCTION 



Superluminal motion observed in Active Galactic Nuclei 
(AGN), especially in the blazars class, seems to be closely 
linked with high-energy emission. Such motion was recently 
observed in the Galaxy (Mirabel & Rodriguez 1994, Hjellm- 
ing & Rupen 1995, Tingay et al. 1995) in the so-called mi- 
croquasars. Nevertheless differences are noticeable in those 
two cases. The latter systems were observed with small value 
of bulk Lorentz factor (around 2.5), while in the former ones 
values of about 10-20 are frequent. It is well known that the 
radiation pressure acting on electron-positron plasma in the 
vicinity of a near Eddington accreting object is very effi- 
cient to accelerate the plasma outwards, since the gravita- 
tional force is around 1000 times weaker than for an electron- 
proton plasma. However, Phinney (1987) has shown that for 
a realistic accretion disc emission, only moderate values of 
bulk Lorentz factors can be reached. Li & Liang (1996), have 
recently proposed that this mechanism could explain the rel- 
atively small 76 ~ 2.5 observed in galactic objects. They 
considered plasma composed with both e + e~ and e~p and 
obtained the equation of motion using the Thomson cross 
section and including gravitation force. To explain higher 
values of 7^, O'Dell (1981) proposed the so-called 'Comp- 
ton rocket' effect, i.e. anisotropic Inverse Compton effect on 
a highly relativistic plasma. It was then argued (Phinney 



1982) that because Compton cooling is always much more 
rapid than bulk acceleration, only small value of 71,00 could 
be reached by this mechanism. However, taking into account 
that in the frame of the 'two-flow' model (Sol, Pelletier & 
Asseo 1989), a pair plasma could be reheated by the tur- 
bulence triggered by a surrounding jet (Henri & Pelletier 
1991), Marcowith, Henri & Pelletier (1995) showed that the 
Compton rocket becomes much more efficient and acceler- 
ates pair plasma to Lorentz factor 75 ~ 10. In another work 
Sikora et al. (1996) studied the radiation drag in AGN jets. 
They included relativistic electron-positron plasma acceler- 
ated in situ and relativistic protons which contribute to the 
inertia of the flow. They considered radiation emitted from 
an accretion disc, partially reprocessed by the outer part 
of the disc or by spherically distributed matter at a given 
distance from the central object. All Compton interactions 
were computed using Thomson cross section. They showed 
that in most cases jets should undergo radiation drag, and 
that the efficiency of this mechanism becomes important for 
purely pair plasma dominating the jet luminosity. The aim 
of this paper is to study how Compton rocket effect can 
accelerate a pure pair plasma in the vicinity of accretion 
disc, taking into account the full KN cross section. Follow- 
ing Marcowith et al. (1995), we consider a blob of e + e~ pair 
plasma with an isotropic energy distribution in the comov- 
ing frame 71^(7') oc 7 /_s , where s is the spectral index. We 
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assume that the acceleration process is efficient enough in 
the jet to get a stationary energy distribution. We assume 
that the dynamics of the pair plasma is decoupled from the 
electron proton component and we do not include gravita- 
tion force. The radiation field is coming from an accretion 
disc. We include KN corrections in the description of the In- 
verse Compton interactions. We study the influence of both 
compactness of the radiation source, spectral index and up- 
per cut-off 'Ymax, and make then comparisons between AGN 
and galactic microquasars. We also consider the influence 
of scattered radiation by a Broad Line Region (BLR) and 
dusty torus around the central black hole. Finally we discuss 
the influence of accretion disc's size that could be relevant 
to the high value of bulk Lorentz factors. 



2 COMPTON ROCKET EFFECT WITH 
KLEIN-NISHINA CORRECTIONS 

2.1 Notations 

All energies are measured in unit m e c 2 . We refer all quanti- 
ties expressed in the blob rest frame by a prime ', all quanti- 
ties in the particle rest frame by a star * and quantities in the 
disc frame are not labelled. Photon energies will be labelled 
by e, and the unit direction vector by fe. We use the KN 
differential cross section (Rybicki & Lightman 1979) given 
by: 



da 



3<7T 



e* + ej 



6 [el 



l+e*(l — cos (f>*) 



ot = ~^-rl is the Thomson cross-section where r e = 
o 

e 2 

is the electron classical radius. This expression 



4neom e c 2 

applies to the scattering of a photon with energy e* and di- 
rection k* in a photon with energy e* and direction fc* , and 
cos <f>* — k*.k*. 



2.2 The general picture 

Figure 1 shows the general configuration of the model. 

The pair plasma is assumed to be described in the bulk 
rest frame by an energy distribution n' e (z,^') oc j'~ s for 



7" 



< 7 < 7" 



with s, 



7" 



and jmax independent of z 



(see paragraph 3.2.1 for a further discussion of this assump- 
tion) . The radiation force is due to soft photons coming from 
a standard accretion disc (Shakura & Sunyaev 1973) around 
a Schwarzschild black hole. The inner radius of the accretion 
disc is n = 3 r g (where r g is the Schwarzschild radius of the 
black hole). The outer radius r e is a free parameter. We use 
the black-body approximation for the disc emission so that 
the specific intensity at a radius r is: 



I„(r)=B„(T eff (r)), 



(1) 



where B v is the Planck function and the effective tempera- 
ture T e f f is given by: 



T ef f(r) 



( 3GMM \ 
\ 8Tvar 3 



1/4 



J [1 - Kn/r) 



1/2] 1/4 



(2) 




Figure 1. The general configuration of the model. Quantities in 
the blob rest frame are denoted by a prime. 



M is the accretion rate, M the mass of the central black 
hole and b is a parameter describing the inner torque on the 
disc (0 < b < 1). In section 4, we also consider the case of a 
two-temperatures disc model (Shapiro, Lightman & Eardley 
1976). This model applies to the hot inner part of the disc 
(r < a few r g ). The specific intensity is then (Sunyaev & 
Titarchuk 1980): 

_ , . exp(— x) /, 1 o 1 3 1 i\ , 

h{r) oc — ^ '- 1 + x + -x 2 + -x + —x 4 - for x > 1 

x \ 2 6 24 / 

I v {r) oc x 4 for x < 1, (3) 

where x = hv/kT e and T e is the electronic temperature given 
by: 



T e (r) = 7.10 8 K 
x{2r/r g ) 1 ' 4 . 



g.s 



1 \ I/ 6 



3M© 



M 



[1 - {n/rf 



/2-1-1/6 



(4) 



In AGN, BLR surrounding the central black hole reprocesses 
a fraction of about 10 per cent of the disc radiation. In para- 
graph 3.2.5, following Sikora et al. (1996), we chose two ex- 
treme cases to describe this emission. The first case is case 
'a' in Sikora et al. (1996) where a fraction \ °f the total disc 
luminosity L d is intercepted by the outer part of the disc 
and reprocessed in Lya lines with a specific intensity given 
by: 



4tt 2 r 



"dr 



ff(f-fLy a ) 



(5) 



The radiation is re-emitted by a ring for n < r < ri. g is a 
rectangular profile with a Doppler width of v/c ~ 2.24 10~ 2 . 
In the second case spherically distributed matter at a dis- 
tance ro from the central engine reemits a fraction \ of the 
total disc luminosity. The emissivity of this radiation source 
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at a distance z on the jet axis is: 



=g(v-VL ya ), (6) 



4 ^oVl+(^) 2 (cos^-l)' 

For z < ro, — 1 < cosd < z/y^z 2 +r% and for z > r , 

\J z 1 — Tq/z < cos# < 1. 

We also study the influence of possible scattered ra- 
diation from hot dust which may be responsible for near 
infrared emission in AGN. The emissivity of such a source 
modeled by a ring located between r$ and ta can be ex- 
pressed as follows : 



x' L d 



Ana J r r4 T 4 (r)rdr 



BA(T(r)), 



(7) 



T is the dust temperature, a is the Stefan-Boltzmann con- 
stant, and x' is tne fraction of the disc luminosity that is 
reprocessed. 



2.3 Computation of the radiative force 



Let 



dp' 



— be the rates of momentum change for one electron 
dt' 

due to Inverse Compton process and /3'c — yjl — 1 /7 /2 c its 
velocity. The force exerced upon the pair plasma is then in 
the blob rest frame: 



(8) 



where u z defines the jet axis direction. This force is accom- 
panied of course by an energy loss that tends to cool the 
plasma. However, as mentioned above, we assume that there 
exists a continuous acceleration mechanism that compensate 
for these losses, so that the particle distribution remains con- 
stant. We assume that this acceleration acts isotropically in 
the plasma frame, and so it does not contribute to the ex- 
pression given in equation (8). We also neglect gravitation 
force, since for near-Eddington objects and cold pairs we 

have — 



^ ~ 10 3 . 

m e 



To calculate more easily F' z , one can use two differ- 
ent approximations: for low energy particles, the KN cross- 
section tends to the Thomson limit and we can apply the 
result of Phinney (1982) for an isotropic pair plasma: 



p/z _ 



(JT 



4nH' 



1 + 3^' 2 ^' 2 ) n ' e ^') d '">''' 



(9) 



where H' is the second Eddington parameter of the radiation 
field (Marcowith et al. 1995, see also Appendix A) defined in 
the blob rest frame. It is linked to the Eddington parameters 
defined in observer's rest frame by: 



H' = 1 ?[(l + $)H-p b (J + K)], 



(10) 



where ftbC is the bulk velocity and 75 is the corresponding 
Lorentz factor. For high energy particles, KN corrections 
become important but we can assume the head-on approxi- 
mation (Blumenthal & Gould 1970) to be valid. In this ap- 
proximation, the velocity of all photons in the particle frame 
is assumed to be anti-parallel to the particle velocity in the 
observer frame. So the momentum loss rate is necessarily in 



the direction of the particle motion. Let x denote this direc- 
tion, the electron energy and momentum loss rate in its rest 
frame are in unit m e c 2 (* refers to this frame): 



dpi 
dt* 

and 



dn* 



da 



dentil 



(e* — el cos <f>*) de*dfli, 



(11) 



(12) 



-IF = - C J dn d^dnt {£ ^ £)d£ ^- 

dn* is the differential photon density in the electron rest 
frame. The momentum change rate in the blob rest frame is 
then given by: 



dpi _ dpi (3' dE* 
dt' dt* c dt* ' 



(13) 



Note that using the KN cross-section, the rate of momentum 
change in the electron rest frame is not zero, contrary to the 
case of Thomson limit or synchrotron emission. One must 
take it into account in computing the energy loss rate in the 
blob rest frame which is not a Lorentz invariant anymore 
(Blumenthal & Gould 1970). In the head-on approximation 
the integral (11) and (12) over e* and SI* can be readily 
performed leading to: 



(14) 



^ = -^t J dn* (/ p ( e *) +/3'f E (e*)). 
The two functions f p and Je are given by: 

/p(0 = (i + 0/(0. 
MO = -e* /(O, 

where the function f(s*) corresponds to the ultrarelativistic 
case (/?' = 1) and is given by: 



/(O = 



1 18e* + 102e* 2 + 186e* 3 + 102e* 4 -20e* 5 
8 e* 3 (l + 2e*) 3 

+ln(l + 2e*)x 
1 (-9 - 141 e* 2 - 60s* - 126e* 3 - 12e* 4 + 24e* 5 ) 
8 e* 3 (l + 2e*) 3 ' 

We then use these expressions in equation (8). 

We can estimate the errors in the two extreme regimes 
described above. In the head-on approximation the first cor- 
rections are roughly — , while in the Thomson limit these 

7 

corrections are ~ 7(e)- Here (e) oc M 1/,4 M _1//2 is the aver- 
age photon energy emitted from accretion disc. We connect 
the two regimes by defining a critical Lorentz factor "i' crit 
for which errors in the head-on approximation (7' > 7 cr ,t) 
and in the Thomson limit (7' < 7 cr it) are of same order. 
This gives j cr it ~ (e)- 1 / 3 . So we may have errors of order 

of 



~ 7crii(e) ~ (e) 2/3 . For AGN (e) - 10~ 4 and we 

find 7crit ~ 20 — 30 with a maximum error ~ 0.2% while for 
a microquasar (e) ~ 10 _ and 7 crit ~ 5 with a maximum 
error ~ 5% . 



-v 2 

I crit 



2.4 Equation of motion 

Following Phinney (1982) we determine the acceleration of 
pair plasma by considering the conservation of stress-energy 
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tensor leading to Phinney's (7) and (8) equations in the bulk 
rest frame: 



(15) 
(16) 



|t HP + Ph»- Ph»] = F'° + b F' z , 
^l( P '+ P 'h b f3 b } = F'*+f3 b F'°. 

Combining these two equations and for a reheated relativis- 
tic plasma (with p' = p'/3) one finds the equation of motion 
(with dz' = fibdt' and z — z')\ 



dz 



(17) 



To compute the radiative force by equation (14) we need the 
differential photon distribution in each electron rest frame. 

dn 

For this we use the Lorentz invariant — (see Blumenthal 

e 

& Gould 1970) and the relation between energy in the ac- 
cretion disc frame and the electron rest frame: 



e* = 76(1 - #,Ms)7(l ~ /V) £ - 



(18) 



In the last equation fj, s is the cosinus of the angle between 
photon direction and jet axis in the accretion disc frame 
(see figure 1) and p! = (i' e (i' s + (1 - ^' e ) 1/2 (l - /4) 1/2 cos<#. 
is the cosinus between the electron and photon direction in 
the bulk rest frame. We use the averaged scattered photon 
energy over the azimutal angle <j> e , and take equation (18) 
with // = iJ,' e iJ,' s . Finally we integrate equation (17) between 
z — 10 r g and z = 10 4 r g (galactic case) or z = 10 5 r g (extra- 
galactic case) for different configurations to determine the 
final bulk Lorentz factor. 



3 RESULTS 

3.1 Equilibrium Lorentz factor 

In any axisymmetric radiation field different from a plane 
wave, there exists an equilibrium Lorentz factor 7(, eg for 
which the radiation force vanishes (-F' z (7f>eq) = 0, see figure 
3). For 7 6 > 7 6eq the pair plasma sees much more photons 
coming forward and then decelerates, while if 76 < jteq the 
plasma is pushed by radiation. The qualitative behaviour of 
a general solution of equation (17) is thus the following: 75 
is first set to the value fbeq- This value increases gradually 
until an altitude z — z cr i t radiation force becomes too weak 
to accelerate the plasma. So it follows a ballistic motion at a 
constant Lorentz factor 7600- This mechanism ensures that 
in most cases the initial value of 75 does not influence the 
terminal value, if the blob is injected bellow z cri t. 

The Thomson regime has been well studied by vari- 
ous authors (Phinney 1982, 1987, Marcowith et al. 1995, 
Sikora et al. 1996). In this case the equilibrium Lorentz fac- 
tor is determined by the simple condition H' — 0, that is 
the radiation flux vanishes in the blob rest frame. Analyti- 
cal calculations in the Thomson regime give the behaviour 
of the equilibrium Lorentz factor and the solutions of the 
equation of motion for a standard accretion disc (Shakura 
& Sunyacv 1973). This yields ~/b eq oc z 1 ^ 4 for z < r e and 
~/beq oc z for z > r e (see Appendix A for this case). Figure 3 
shows the function 7(, e9 obtained in the Thomson limit for 
different values of r e in comparison with these two asymp- 
totic regimes. Note that for decreasing size of the accretion 
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Figure 2. The radiative force as a function of the bulk Lorentz 
factor 7;, at different altitudes from the central black hole. 

disc, the radiation field is more and more anisotropic and the 
function jb eq increases more rapidily with z. For z cr it < r E 

f(l' 2 )lY /7 

one finally gets for a hot plasma j boo oc ±j — ■— , while 



for z cr it > r e , 7f>oc oc 



M 3/4 <7 ,2 V X : 



(see Appendix 



4irm e c 3 ri 



is the compactness of the source, and 



UlW x / /n' e ( 7 ')d7' 



A). / 

{/) = 

Figure 4 shows the equilibrium Lorentz factor calcu- 
lated including KN corrections for different configurations. 
We compare the results with the correspondant Thomson 
solution. The effect of the KN corrections is to reduce the 
contribution of high energetic collisions (e* > 1) in the net 
radiation force seen by the plasma. Since the more energetic 
photons come from the inner part of the accretion disc, in 
the blob rest frame they contribute to accelerate the plasma. 
On the other hand the outer and colder part contribute to 
decelerate it. We then expect that the equilibrium Lorentz 
factor is reduce when including KN corrections as shown in 
figure 4. The importance of the difference depends on the 
particle and photon distributions. It is predominant for mi- 
croquasars in which accretion disc radiates more energetic 
photons than in the extragalactic case. 

3.2 Influence of the different parameters 

3.2.1 General considerations 

To determine the terminal Lorentz factor for different con- 
figurations, we solve numerically the equation of motion for 
various radiation fields and particle distributions. We show 
in figure 5 typical solutions in the case of a stellar black 
hole. The terminal value of the bulk Lorentz factor does not 
depend on the initial value of jt init as discussed in para- 
graph 3.1. The general behaviour of the solution discussed 
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Figure 3. Equilibrium Lorentz factor fbeq as a function of z in 
the Thomson limit. lb eq is given for r e = 10 r 9 , r e = 10 2 r g , 
r e = 10 3 r g and r e = 3.10 3 r g . We also represent the two asymp- 
totic regimes jt> eg oc z 1 / 4 for 2 < r e and jt, eq oc z for z > r e , with 
r e = 3.10 3 r o . 
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Figure 4. The Equilibrium Lorentz factor calculated including 
KN corrections compared to Thomson solution. The figure cor- 
responds to a stellar black hole M = 5Mq with y m ax = 10 7 
for spectral index s = 1.5 and s = 5. The external radius is 
r e = 3 10 3 r g and L = L Edd . 

in this paragraph still holds even including KN corrections. 
The critical point z cr it is reached rather close to the central 
engine (before 10 4 Schwarzchild radii) . For z > z cr u the mo- 
tion is nearly balistic and so independent of the radiation 
force which has become too weak. It is so independent of 
any variation of the pair distribution unless these variations 
strengthen the radiation force. This scenario would require 
a more efficient acceleration mechanism when moving away 



from the central source, which is very unlikely. Therefore 
our assumption of a stationary pair energy distribution on 
a large range of z does not influence strongly the terminal 
value of 7b, or in other words this value is essentially deter- 
mined by the local parameters at the critical distance. 

3.2.2 Influence of the energy upper cut-off and spectral 
index 

Figure 6 illustrates the influence of the spectral index (for 
1.5 < s < 5) and the energy cut-off (for 10 3 < -y ma x < 
10 7 ) on the terminal Lorentz factor. We chose M — 5M@ 
representative for stellar black holes and M — 1O 9 M0 for 
supermassive black holes. The calculations were carried out 
for L = LEdd and L = Q.lLsdd, where L is the luminosity 
of the accretion disc. 

Results are very sensitive to the spectral index value. 
There are 3 different behaviours according to the value of s. 

(i) for s < 2 there exists a maximum terminal Lorentz 
factor as a function of y ma x and the curve intercepts the 
other ones. 

(ii) for 2 < s < 3 there still exists a maximum, but 
less pronounced. The variation as a function of y ma x arc 
smoother. 

(iii) for s > 3 no variation with j ma x- 

We find that for low value of 7max, our solutions agree with 
the Thomson regime solutions. Nevertheless KN corrections 
reduce the efficiency of the Compton rocket effect. As a 
matter of fact, in the Thomson regime, an increase of 7max 
leads to an increase of (j 12 ) (for s < 3) and so of 7600. This 
mechanism is valid until KN corrections begin to dominate, 
roughly when j max (e) ~ 1. So when j ma x is greater than 
(e) -1 the radiation force does not increase anymore whereas 
the plasma inertia p' is much more important. This leads to 
a less efficient acceleration mechanism. This effect is larger 
for small indexes explaining the inversion of the curve for 
high 7max- Indeed we find that acceleration is much more 
efficient for s = 2 and j ma x = 10 7 than for s = 1.5 and the 
same value of f m ax, in the case of stellar black hole. When 
steepening the pair distribution, the radiation force is dom- 
inated by the low energy part of the distribution. This fact 
explains why no variation is apparent with jmax for s > 3. 
The plasma behaves dynamically as a cold one, and we find 
small value of terminal Lorentz factor. Finally including KN 
corrections in the calculations gives rise to an absolute upper 
limit to maximal Lorentz factor for a given luminosity. 

3.2.3 Influence of the black hole mass 

As discussed above the influence of the mass of the central 
black hole is predominant. Stellar black holes with soft X-ray 
emission ((e) ~ 1CP 2 ) are less efficient in accelerating blob 
of pair plasma than supermassive black hole with softer ra- 
diated emission ((e) ~ 1CF 4 ), because KN saturation effects 
occur at much lower energy. A more realistic description of 
the accretion disc around stellar black holes reinforces this 
discrepency. As shown in figure 7 the radiation emitted from 
a two-temperature disc (Shapiro et al. 1976) leads to smaller 
7i,oo than in the case of standard accretion disc radiation. 
Because it is well established that accretion disc around stel- 
lar black holes should radiate up to a few keV (as in a two- 
temperatures disc), KN corrections play an important role 
in this case. 
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Figure 5. Solutions of the equation of motion in the case of a 
stellar black hole (M = 5Mq). We chose two different initial con- 
ditions 7 6 init = 2 and 7 fc init = 5. r e = 3 10 3 r g and L = L Edd . 



3.2.4 Influence of the luminosity 

As (e) tx M 1 / 4 M~ 1/ ' 2 , the luminosity of the disc also directly 
infiuences the maximum of the function 7(,oo as a function 
of jmax- As shown in figure 6, the maximum takes place 
at higher j max when the luminosity decreases. Besides less 
luminous systems contribute to a lower radiation force and 
so to a less efficient acceleration. In the Thomson regime one 
has a dependance 7^ oc L 1 ^ 7 . 

3.2.5 Effect of scattered radiation 

All the results described above are obtained when studying 
the disc radiation alone. We also include in our calculation 
BLR radiation fields corresponding to two cases: 

a) re-emission from a ring located between n and r2 and 
with an emissivity given by equation (5) 

b) re-emission from spherically distributed matter at a 
distance ro from the central black hole with an emissivity 
given by equation (6). 

Figure 8 displays the equilibrium Lorentz factor in the pres- 
ence of a BLR located between n = 1 4 r g and ri = 1 5 r g 
(case a), figure 8) in the Thomson regime. We also plotted 
the equlibrium Lorentz factor including KN corrections as 
well as the solution of the equation of motion for a plasma 
with s = 2 and jmax = 10 5 . As one can see, the effect of 
BLR on 76 e q is very weakened by KN corrections. This can 
be understood because the photons coming from the BLR 
are blueshifted by the relativistic motion in the blob rest 
frame whereas the photons coming from the disc are red- 
shifted. So the dragging force from the BLR is much more 
reduced by KN corrections than the accelerating one due 
to disc photons. Due to the weakness of the radiation field, 
the dynamical solution j b (z) is still less affected than the 



equilibrium value. The case of a spherical shell located at 
ro = 1 4 r g (case b)) is illustrated in figure 9. One can see 
that the diffused radiation field strongly affects the motion 
that is almost stopped at the crossing of the shell, where ra- 
diation density is dominated by the isotropically scattered 
photons. However, the plasma is quickly reaccelerated after 
the crossing, and high Lorentz factor can be reached again. 

Sikora et al. (1996) argued that scattered radiation 
should lead to efficient radiation drag. However, the influ- 
ence of scattered radiation is strongly governed by the posi- 
tion of scattering clouds with respect to the critical distance 
z cr it : if they lie below this critical distance, the plasma 
will be temporarily braked during the crossing of scatter- 
ing region, but will be quickly reaccelerated after it. If the 
clouds lie above z cr it, the terminal Lorentz factor can indeed 
be strongly affected. Figure 10 and II display the terminal 
Lorentz factor as a function of ^ m ax for different distances 
of the BLR. It can be seen that if it is close enough, the BLR 
can give still higher Lorentz factor than for the disc alone 
for the highest value of 7 max . This is because the solid angle 
substended by it at z crit is so small that its radiation field 
has an accelerating rather than decelerating effect. But even 
for a BLR between 10 4 and 10 s r g , terminal Lorentz factor 
around 10 are clearly reachable. We conclude that the pres- 
ence of broad lines can affect 7600 , but does not prevent in 
general highly relativistic motions. 

We also discuss the effect of a possible emission from a 
hot dusty torus surrounding the central black hole outside 
the BLR region. The effect of such a scattered radiation 
field is much more important because KN corrections are 
here almost negligible. Figure 12 shows the terminal Lorentz 
factor in the presence of an IR emitting ring heated at T = 
1500AT, located between r$ and r±. It is obvious that the 
terminal Lorentz factor is strongly reduced to values around 
3 for the less favourable case (r 3 = 10 4 r 9 , r 4 = 10 5 r 9 ). The 
situation is a little less dramatic for more distant sources 
(rz = 10 5 r 9 , T4, — 10 6 r g ) because the isotropic radiation 
density is lowered, and 71,00 of 5 can be reached. 

3.2.6 Influence of the accretion disc size 

For small value of r e the radiation is more anisotropic 
and so more efficient to accelerate the pair plasma (see figure 
13). In this figure r e = 10 r g and we show the solutions for 
which we obtain the largest 71,00 in extragalactic and galac- 
tic cases. In this configuration the radiation force is more 
efficient and 7600 can be as high as 60 in the extragalac- 
tic case. The dependance on s and j max is shown in figure 
14, where we extend the previous calculation to the case 
r e = 10 r g . We find the same global behaviour of 7600 with 
spectral index and y ma x. 



4 COMPARISONS WITH OBSERVATIONS 

Our model gives good agreement with observations of rela- 
tivistic blob ejection in AGN and microquasars for the 'disc 
alone' solution. Figure 6 shows that in the most favourable 
configuration (s = 1.5, "y ma x ~ 10 6 and L = L Edd ) Compton 
rocket effect is able to accelerate pair plasma up to Lorentz 
factor 71,00 > 20. This value is only gradually reached, the jet 
being much slower at small distances. One can obtain higher 
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Figure 6. Terminal Lorentz factor -ytoa as a function of f ma x for different value of spectral index s. The two top panncls correspond to 
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value for super-Eddington systems. One can notice that ob- & Cohen 1994). 

servations of the faster extragalactic superluminal motion This value is however strongly dependent on the spec- 

correspond to such value of bulk Lorentz factor (Vermeulen tal index s and the high energy particle distribution cut-off 
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Figure 8. Equilibrium Lorentz factor as a function of z for radia- 
tion coming from an accretion disc and BLR radiation from a ring 
located between n = 10 4 r 9 and T2 = 10 5 r 9 . The emissivity is de- 
scribed by equation (5) with \ = 0.1 and a = 2. We also show 
the solution of the equation of motion for this case. The solid 
curve shows "fteq m Thomson regime. The plasma parameters are 
s = 2 and ^ m ax = 10 5 . The black hole mass is M = 1O 9 M and 
L = L Edd . 



1 1 

7b a, M=10 9 M a s=2 7moI =10 5 





i - 


— - A M='°\> s=2 7mo> =10 5 




7 

.7 


7b eq Thomson regime 




/ - 


case b} 






i i 


\ if 





10 

z / r „ 



Figure 9. Same as figure 8 for a BLR radiation from a shell 
located at rrj = 10 4 r g . The emissivity is described by equation 
(6) with x = 0.1. 

7 mol . The precise value of s is not obvious to derive from 
observations. High energy spectra show typical X-ray spec- 
tral indexes around 0.5, which correspond to s ~ 2. There is 
very often a spectral break in the MeV range: this could be 
attributed to a break in the particle distribution, that would 
correspond to jmax = W 2 . However, Marcowith et al. (1995) 
have shown that this break could be very well reproduced 



Figure 10. Terminal Lorentz factor "ftoo as a function of ^/ ma x 
for the spectral index s = 2. The solid line shows the solution 
obtained including KN corrections for a standard accretion disc. 
The scattered radiation is described by equation (5) with \ = 0.1, 
a = 2 and different locations of the BLR ring. 
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Figure 11. Same as figure 10 for scattered radiation described 
by equation (6) with x = 0.1 and different locations of the BLR 
shell. 

by 7 — 7 absorption with an actual particle distribution giv- 
ing the primary photon spectrum that can extend to much 
higher energy. So high energy spectra may not be good indi- 
cators of the upper cut-off y ma x ■ Moreover, the high energy 
emission probably takes place at relatively small distances 
(around 10 2 r 9 ) from the center, well below z cri t: the final 
bulk Lorentz factor is not yet reached at this distance. On 
the other hand, the detection of photons with at least 30 
GeV, and even above TeV for some BL Lacs, implies an 
upper cut-off j max > 10 5 . 
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Figure 12. The influence of a dusty ring combined with a BLR 
on the terminal Lorentz factor 7^00 as a function of -f m ax , for the 
spectral index s = 2. The solid line shows the solution obtained 
including KN corrections for a standard accretion disc. The BLR 
ring emissivity is described by equation (5) with x = 0.1, a = 2. 
The dusty ring emissivity is given by equation (7) with \' = 0.1 or 
x' = 0.05. We chose different locations for these two components. 



- 7 b eq Thomson regime / 






// 

/'/ y b M=10 \ s= 


'■5 7™=10 6 -- 


--- y b init-4 


\ /" 






1 / 7b M =5M s= 


■5 ? m „„=10 4 


- - y b init-4 



Figure 13. Solutions of the equation of motion for a supcrmassivc 
black hole (M = 1O 9 M ) and a stellar black hole (M = 5M Q ) 
with a compact accretion disc. r e = 10 r g and L = L^dd- lb * n ** 
is the initial condition of the solution. We also represent the equi- 
librium Lorentz factor in Thomson regime . 

Radio spectral indexes are also difficult to assess be- 
cause of synchrotron self-absorption, especially for radio- 
fiat quasars where the observed spectrum results probably 
from the superposition of many self-absorbed spectra. The 
optically thin index ranges mostly from 0.5 to 1, which 
corresponds to 2 < s < 3. With reasonable parameters 



(2 < s < 3 and ~{ m ax < 10 5 ), our model predicts typical 
values 4 < 71,00 < 10 which are the most frequently ob- 
served (Vermeulen & Cohen 1994). Moreover, as shown in 
Sect. 3.2.5, KN effects prevents from strong Compton drag 
induced by BLR radiation in the vicinity of the central black 
hole. The fastest superluminal motions may be attributed to 
those objects for which the BLR is weak and/or closest to 
the central object. 

On the other hand, the slowest motions (7600 < 4) can 
be obtained in the presence of a rich and extended environ- 
ment of scattering material, such as BLR clouds and dusty 
torus. The presence of dust is inferred from an enhance- 
ment in the near infrared in some quasars spectra (Bairvai- 
nis 1987). Nevertheless a nonthermal infrared component, 
attributed to synchrotron radiation from the relativistic jet, 
is also usually observed in radio- loud AGN. This compo- 
nent is generally predominant in flat spectrum radio quasars 
(Neugebauer et al. 1986) and is necessary to explain rapid 
variations observed in infrared flaring objects. We can specu- 
late that these objects with the largest superluminal motion 
are those where the scattered thermal component is par- 
ticularly weak. Such 'nonthermal' objects should not suffer 
strong radiation drag, while 'thermal' ones should have the 
lowest terminal Lorentz factors. In conclusion the diversity 
among observed superluminal motions can be reproduced 
by our model by considering both influence of the plasma 
acceleration mechanism and AGN environment. 

The observed value of Lorentz factor of about 2.5 for 
the two microquasars (GRS1915+105, Mirabel & Rodriguez 
1994 and GROJ1655-40 , Hjemlling & Rupen 1995) with 
large spectral indexes (respectively s ~ 4, Finoguenov et 
al. 1994 and s ~ 4.6, Harmon et al. 1995) are consistent 
with our results. We show figure 15 the dependance of 71,00 
on compactness of the source for s — 4 and s — 5. We 
find that observations require L ~ 0.2 — 0.3L E dd, which is 
very close to the result by Li & Liang (1996). This result 
is not strongly affected by the mass of the central black 
hole if the compactness of the source remains the same. The 
steep spectrums observed in these two objects ensure us that 
the terminal Lorentz factor is just dependant on the low 
energy cut-off of the electron distribution. In this case a more 
realistic description of the accretion disc emission does not 
change our results for such spectral indexes (figure 7). Yet 
the maximal value of 7^00 is model dependant. The part of 
high energy emission contributes to decrease this value from 
the one obtained for a standard accretion disc. We find that 
value of about 6 can be reached in the most favourable case. 

Finally, 7600 is much higher for a small sized accretion 
disc. We find value of about 60 in the case of an external 
radius of 10 r g in the case of supermassive black hole. 



5 CONCLUSION 

We have considered the bulk acceleration of an electron- 
positron pair plasma in the vicinity of a central black hole. 
The acceleration is due to Compton rocket effect on the 
plasma and the radiation force originates from a standard 
accretion disc emission. The pair plasma is continuously re- 
heated by an efficient turbulent mechanism wich takes place 
in the frame of the 'two-flow' model. Thus we assumed in 
our calculations a stationary power-law energy distribution 
for the pair. We included KN corrections in the computation 
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Figure 14. Terminal Lorentz factor 7^ as a function of -f max for different value of spectral index s. The two pannels correspond to 
r e = Wr g , L = L Edd , M = lO 9 M (left) and M = 5M Q (right). 
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Figure 15. The terminal Lorentz factor as a function of L/LEdd 
and for two spectral indexes (s = 4 and s = 5). M = 5Mq, 
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of the radiation force, and solved numerically the equation 
of motion. We studied configurations relevant to relativistic 
motion in AGN and galactic microquasars. The main results 
of our calculations can be summarized as follows: 

(i) For a given luminosity, the terminal Lorentz factor 7^ 
admits an absolute maximum due to KN corrections. Val- 
ues of about 20 can be reached in the extragalactic case, for 
sufficiently flat spectrum (s ~ 1.5) and accretion at the Ed- 
dington rate which may correspond to the highest observed 
relativistic motion. For more reasonable plasma parameters 
(s < 2 and ^ ma x < 10 6 ) the Compton rocket effect can ac- 
count for the typical value of terminal Lorentz factor inferred 
from observations (7600). 



(ii) Scattered radiation from extended BLR or dusty torus 
can brake efficiently relativistic motion, even for a high en- 
ergy plasma. This Compton drag and a weak plasma heat- 
ing can be responsible for the lowest terminal Lorentz fac- 
tor observed. The highest superluminal motion could be at- 
tributed to objects with particularly weak diffuse component 
and very efficient heating. 

(iii) For stellar black holes KN corrections are important 
leading to smaller values than for a supermassive black hole. 
Recent observations of relativistic ejection in galactic micro- 
quasar are consistent with our results. 
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APPENDIX A: EQUILIBRIUM LORENTZ 
FACTOR AND TERMINAL LORENTZ FACTOR 
IN THOMSON REGIME 

Eddington parameters are defined as follows: 
J = 4^ y hiP)AQ,Au, 



K=^- J u 2 I„(n)dQdv. 



(Al) 



Marcowith et al. (1995) show that these parameters can be 
expressed in the case of standard accretion disc (Shakura & 
Sunyaev 1973) as a function of the integral: 

I{a,M) = j e [l/« a (l + £ V) /3 ][l-(l/ W ) 1/2 ]^d W; (A2) 

with u — r/ri, u e — r e /ri and e = Vi/z. They obtain: 



J = (L/87v 2 rf)(l - 2&/3)- 1 e 2 /(2,3/2, 1), 
H = (L/8ir 2 r?)(l - 2fc/3)- V/(2, 2, 1), 
K = (L/8n 2 r 2 )(l - 2&/3) _1 e 2 J(2, 5/2, 1). 



(A3) 



They study the case n < z < r e for which s < 1 and u > 1. 
They find that the equilibrium Lorentz factor is given by: 



"/beg ~ £ 



1/4 _ Z 



1/4 
,V4 



(A4) 



In the case z > r e we have eu = r /z < 1 for 1 < u < u e . So 
we can expand the term (1 + e 2 u 2 )~^ into series: 

(1 + e 2 u 2 )^ = ^(~l) fc ^ ± V-jf + k ' 1] £ 2 V fc . (A5) 

k 

Hence in the case a = 2, 7 = 1 and integrating over u ,we 
obtain: 



J(2,/3,l)=5> 2fc (-l 



fc /3Q3 + l)...Q3 + fc~l) 
fc! 



u?*- 1 - 1 u? fe - 3/2 - 1 



2fe - 1 2fc - 3/2 

The first coefficients in e are then: 



1(2,(1,1) = A ~ f3A l£ 2 + 



with 



f3([] + i) 



A 2 s\ 



(A6) 



(A7) 



, 1 -1,2 -3/2 1 
3 ~ Ue 3 Ue ~ 3' 

Al — U e — 2u~J 2 + 1 ~ W e , 

tt 3 2lte /2 1 U\ 

A 2 = — 1 ~ — . 

3 5 15 3 



(A8) 



The equivalents are given for u e 3> 1. This gives for the 
Eddington parameters: 

J = (L/87r 2 r 2 )(l - 2b/3)- 1 e 2 (A - 3/2A 1 e 2 + 15/8A 2 e 4 ), 
H = (L/87r 2 r 2 )(l - 2b/zy 1 e 2 {A - 2A x e 2 + 3A 2 e 4 ), (A9) 
K = (L/87r 2 r 2 )(l - 2&/3) _1 e 2 (A - 5/2^ie 2 + 35/8A 2 £ 4 ). 
The equilibrium Lorentz factor is 7 [, eq = (1— /3 2 eq )~ 1//2 where 



x — (x — 1) 



1/2 



with 



J + K 

x = ^2iT- 

Using equations (A9) we find with u e 3> 1: 



1 1 3 4 

x = 1 + -u e e , 



Jbeq 



3/2 
U e 2 



and finally 

-3/4 
-fbeq =U e £ — 



3/4 1/4 ' 

r e r- 



(A10) 
(AH) 

(A12) 
(A13) 

(All) 



The asymptotic solution of equation (17) is approxima- 
tively 7600 = ~y b eq(z cr it), where z cr it is the location where 
the radiative force becomes to weak to maintain "f b (z) = 
j beq (z). This occurs when the evolution of j beq , i.e. Azo = 
Jbeq/id-ybeq/dz) = dlog(«) / (zdlog^beq ) ) , is larger than the 
evolution of -y b towards j beq , i.e. A21 = (dy b /dz) / (y b - y beq ) . 
Using equation (17) one finds: 
\ -1 



/ 1 dF' z 



(A15) 



In the Thomson regime we can use equation (9) and (10) 
which give in the relativistic case: 



dF' z a T 8tt. , 2 .dH' 2H n 

— = 5- (7 H — = — (l-x/Pb) 

d-y b c 3 d 7fc 7 6 



(A16) 



With (A9), (A12), (A13), (A14) and ignoring some terms of 
order unity one finally gets: 

y' 2 )V /4 

■■<■• < - I ' - ' , r-T^rr I , (A 17) 



m e c 3 (7') 



and, 

Iboo ' 



3/4 (Y 2 ) 

(70 



1/4 



(Ais) 



In the case n < z < r e , using equation (A4) one finds 
La T <7' 2 > 



Zcrit ~ T% 



and finally, 



m e c 3 r 4 (7') 



76o 



<7' ; 



<7'> 
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(A19) 
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